The macroevolutionary singularity of snakes.

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

FieldPrism: A system for creating snapshot vouchers from field images using photogrammetric markers and QR codes.

Premise: Field images are important sources of information for research in the natural sciences. However, images that lack photogrammetric scale bars, including most iNaturalist observations, cannot yield accurate trait measurements. We introduce FieldPrism, a novel system of photogrammetric markers, QR codes, and software to automate the curation of snapshot vouchers. Methods and Results: Our photogrammetric background templates (FieldSheets) increase the utility of field images by providing machine-readable scale bars and photogrammetric reference points to automatically correct image distortion and calculate a pixel-to-metric conversion ratio. Users can generate a QR code flipbook derived from a specimen identifier naming hierarchy, enabling machine-readable specimen identification for automatic file renaming. We also developed FieldStation, a Raspberry Pi-based mobile imaging apparatus that records images, GPS location, and metadata redundantly on up to four USB storage devices and can be monitored and controlled from any Wi-Fi connected device. Conclusions: FieldPrism is a flexible software tool designed to standardize and improve the utility of images captured in the field. When paired with the optional FieldStation, researchers can create a self-contained mobile imaging apparatus for quantitative trait data collection.

From leaves to labels: Building modular machine learning networks for rapid herbarium specimen analysis with LeafMachine2.

Premise: Quantitative plant traits play a crucial role in biological research. However, traditional methods for measuring plant morphology are time consuming and have limited scalability. We present LeafMachine2, a suite of modular machine learning and computer vision tools that can automatically extract a base set of leaf traits from digital plant data sets. Methods: LeafMachine2 was trained on 494,766 manually prepared annotations from 5648 herbarium images obtained from 288 institutions and representing 2663 species; it employs a set of plant component detection and segmentation algorithms to isolate individual leaves, petioles, fruits, flowers, wood samples, buds, and roots. Our landmarking network automatically identifies and measures nine pseudo-landmarks that occur on most broadleaf taxa. Text labels and barcodes are automatically identified by an archival component detector and are prepared for optical character recognition methods or natural language processing algorithms. Results: LeafMachine2 can extract trait data from at least 245 angiosperm families and calculate pixel-to-metric conversion factors for 26 commonly used ruler types. Discussion: LeafMachine2 is a highly efficient tool for generating large quantities of plant trait data, even from occluded or overlapping leaves, field images, and non-archival data sets. Our project, along with similar initiatives, has made significant progress in removing the bottleneck in plant trait data acquisition from herbarium specimens and shifted the focus toward the crucial task of data revision and quality control.

Fish Welfare in Public Aquariums and Zoological Collections.

A wide variety of fish species have been displayed in public aquariums and zoological collections for over 150 years. Though the issue of pain perception in fish is still being debated, there is no disagreement that negative impacts on their welfare can significantly affect their health and wellbeing. A general description of the basic biological requirements for maintaining fish in captive environments is presented, but species-specific information and guidelines should be developed for the multitude of species being maintained. A combination of behavioral, performance, and physiological indicators can be used to assess the well-being of these animals. Ultimately, the goal for optimizing the welfare of fish should be to provide the best possible environment, husbandry, and social interactions to promote natural species-specific behaviors of the fish in captivity.

Global hotspots of plant phylogenetic diversity.

Regions harbouring high unique phylogenetic diversity (PD) are priority targets for conservation. Here, we analyse the global distribution of plant PD, which remains poorly understood despite plants being the foundation of most terrestrial habitats and key to human livelihoods. Capitalising on a recently completed, comprehensive global checklist of vascular plants, we identify hotspots of unique plant PD and test three hypotheses: (1) PD is more evenly distributed than species diversity; (2) areas of highest PD (often called 'hotspots') do not maximise cumulative PD; and (3) many biomes are needed to maximise cumulative PD. Our results support all three hypotheses: more than twice as many regions are required to cover 50% of global plant PD compared to 50% of species; regions that maximise cumulative PD substantially differ from the regions with outstanding individual PD; and while (sub-)tropical moist forest regions dominate across PD hotspots, other forest types and open biomes are also essential. Safeguarding PD in the Anthropocene (including the protection of some comparatively species-poor areas) is a global, increasingly recognised responsibility. Having highlighted countries with outstanding unique plant PD, further analyses are now required to fully understand the global distribution of plant PD and associated conservation imperatives across spatial scales.

Compositional shifts associated with major evolutionary transitions in plants.

Heterogeneity in gene trees, morphological characters, and composition has been associated with several major plant clades. Here, we examine heterogeneity in composition across a large transcriptomic dataset of plants to better understand whether locations of shifts in composition are shared across gene regions and whether directions of shifts within clades are shared across gene regions. We estimate mixed models of composition for both nucleotide and amino acids across a recent large-scale transcriptomic dataset for plants. We find shifts in composition across both nucleotide and amino acid datasets, with more shifts detected in nucleotides. We find that Chlorophytes and lineages within experience the most shifts. However, many shifts occur at the origins of land, vascular, and seed plants. While genes in these clades do not typically share the same composition, they tend to shift in the same direction. We discuss potential causes of these patterns. Compositional heterogeneity has been highlighted as a potential problem for phylogenetic analysis, but the variation presented here highlights the need to further investigate these patterns for the signal of biological processes.

Identifying Climatic Drivers of Hybridization with a New Ancestral Niche Reconstruction Method.

Applications of molecular phylogenetic approaches have uncovered evidence of hybridization across numerous clades of life, yet the environmental factors responsible for driving opportunities for hybridization remain obscure. Verbal models implicating geographic range shifts that brought species together during the Pleistocene have often been invoked, but quantitative tests using paleoclimatic data are needed to validate these models. Here, we produce a phylogeny for Heuchereae, a clade of 15 genera and 83 species in Saxifragaceae, with complete sampling of recognized species, using 277 nuclear loci and nearly complete chloroplast genomes. We then employ an improved framework with a coalescent simulation approach to test and confirm previous hybridization hypotheses and identify one new intergeneric hybridization event. Focusing on the North American distribution of Heuchereae, we introduce and implement a newly developed approach to reconstruct potential past distributions for ancestral lineages across all species in the clade and across a paleoclimatic record extending from the late Pliocene. Time calibration based on both nuclear and chloroplast trees recovers a mid- to late-Pleistocene date for most inferred hybridization events, a timeframe concomitant with repeated geographic range restriction into overlapping refugia. Our results indicate an important role for past episodes of climate change, and the contrasting responses of species with differing ecological strategies, in generating novel patterns of range contact among plant communities and therefore new opportunities for hybridization. The new ancestral niche method flexibly models the shape of niche while incorporating diverse sources of uncertainty and will be an important addition to the current comparative methods toolkit. [Ancestral niche reconstruction; hybridization; paleoclimate; pleistocene.].

The phylogeny and global biogeography of Primulaceae based on high-throughput DNA sequence data.

The angiosperm family Primulaceae is morphologically diverse and distributed nearly worldwide. However, phylogenetic uncertainty has obstructed the identification of major morphological and biogeographic transitions within the clade. We used target capture sequencing with the Angiosperms353 probes, taxon-sampling encompassing nearly all genera of the family, tree-based sequence curation, and multiple phylogenetic approaches to investigate the major clades of Primulaceae and their relationship to other Ericales. We generated dated phylogenetic trees and conducted broad-scale biogeographic analyses as well as stochastic character mapping of growth habit. We show that Ardisia, a pantropical genus and the largest in the family, is not monophyletic, with at least 19 smaller genera nested within it. Neotropical members of Ardisia and several smaller genera form a clade, an ancestor of which arrived in the Neotropics and began diversifying about 20 Ma. This Neotropical clade is most closely related to Elingamita and Tapeinosperma, which are most diverse on islands of the Pacific. Both Androsace and Primula are non-monophyletic by the inclusion of smaller genera. Ancestral state reconstructions revealed that there have either been parallel transitions to an herbaceous habit in Primuloideae, Samolus, and at least three lineages of Myrsinoideae, or a common ancestor of nearly all Primulaceae was herbaceous. Our results provide a robust estimate of phylogenetic relationships across Primulaceae and show that a revised classification of Myrsinoideae and several other clades within the family is necessary to render all genera monophyletic.

Analysis of the core bacterial community associated with consumer-ready Eastern oysters (Crassostrea virginica).

Shellfish, such as the Eastern oyster (Crassostrea virginica), are an important agricultural commodity. Previous research has demonstrated the importance of the native microbiome of oysters against exogenous challenges by non-native pathogens. However, the taxonomic makeup of the oyster microbiome and the impact of environmental factors on it are understudied. Research was conducted quarterly over a calendar year (February 2020 through February 2021) to analyze the taxonomic diversity of bacteria present within the microbiome of consumer-ready-to-eat live Eastern oysters. It was hypothesized that a core group of bacterial species would be present in the microbiome regardless of external factors such as the water temperature at the time of harvest or post-harvesting processing. At each time point, 18 Chesapeake Bay (eastern United States) watershed aquacultured oysters were acquired from a local grocery store, genomic DNA was extracted from the homogenized whole oyster tissues, and the bacterial 16S rRNA gene hypervariable V4 region was PCR-amplified using barcoded primers prior to sequencing via Illumina MiSeq and bioinformatic analysis of the data. A core group of bacteria were identified to be consistently associated with the Eastern oyster, including members of the phyla Firmicutes and Spirochaetota, represented by the families Mycoplasmataceae and Spirochaetaceae, respectively. The phyla Cyanobacterota and Campliobacterota became more predominant in relation to warmer or colder water column temperature, respectively, at the time of oyster harvest.

Lagrange-NG: The next generation of Lagrange.

Computing ancestral ranges via the Dispersion Extinction and Cladogensis (DEC) model of biogeography is characterized by an exponential number of states relative to the number of regions considered. This is because the DEC model requires computing a large matrix exponential, which typically accounts for up to 80% of overall runtime. Therefore, the kinds of biogeographical analyses that can be conducted under the DEC model are limited by the number of regions under consideration. In this work, we present a completely redesigned efficient version of the popular tool Lagrange which is up to 49 times faster with multithreading enabled, and is also 26 times faster when using only one thread. We call this new version Lagrange-NG (Lagrange-Next Generation). The increased computational efficiency allows Lagrange-NG to analyze datasets with a large number of regions in a reasonable amount of time, up to 12 regions in approximately 18 min. We achieve these speedups using a relatively new method of computing the matrix exponential based on Krylov subspaces. In order to validate the correctness of Lagrange-NG, we also introduce a novel metric on range distributions for trees so that researchers can assess the difference between any two range inferences. Finally, Lagrange-NG exhibits substantially higher adherence to coding quality standards. It improves a respective software quality indicator as implemented in the SoftWipe tool from average (5.5; Lagrange) to high (7.8; Lagrange-NG). Lagrange-NG is freely available under GPL2. [Biogeography; Phylogenetics; DEC Model.].

First isolation of Carnobacterium maltaromaticum from farmed Rainbow Trout in Virginia.

OBJECTIVE: Carnobacterium maltaromaticum is considered an emerging pathogen of salmonids in the United States and around the world. METHODS: Bacterial cultures obtained from the posterior kidney and skin of moribund Rainbow Trout Oncorhynchus mykiss from a commercial aquaculture facility in Virginia, USA, grew C. maltaromaticum, which was confirmed by additional phenotypic and molecular characterization. RESULT: A presumptive diagnosis based on the clinical signs, necropsy observations, histopathology, and bacterial cultures was bacterial septicemia due to C. maltaromaticum. CONCLUSION: This represents the first documentation of C. maltaromaticum in Rainbow Trout from Virginia.

Development of a Controlled Laboratory-scale Inoculation System to Study Vibrio parahaemolyticus-oyster Interactions.

Prevalence of seafood-borne gastroenteritis caused by the human pathogen Vibrio parahaemolyticus is increasing globally despite current preventative measures. The United States Centers for Disease Control have designated V. parahaemolyticus as a reportable emerging human pathogen. The Eastern oyster (Crassostrea virginica) is a natural reservoir of the bacterium in marine environments, but little is actually known regarding interactions between oysters and V. parahaemolyticus. Therefore, a laboratory-scale Biosafety Level-2 (BSL2) inoculation system was developed wherein Chesapeake Bay region oysters harvested during summer or winter months, were exposed to the clinical RIMD2210633 strain carrying a chloramphenicol-selective marker (VP RIMDmC). Homogenized whole oyster tissues were spread on selective and differential agar medium to measure viable VP RIMDmC levels. Endogenous Vibrio spp. cell numbers were significantly reduced followed chloramphenicol treatment and this likely contributed to higher VP RIMDmC oyster-associated levels, especially using winter-harvested animals. Summer-harvested oysters had significantly higher existing Vibrio levels and a lower level of artificial oyster-associated VP RIMDmC. Thus, the pre-existing microbiome appears to afford some protection from an external V. parahaemolyticus challenge. Overall, this system successfully enabled controlled manipulation of parameters influencing V. parahaemolyticus-oyster interactions and will be useful in safely testing additional pertinent environmental variables and potential mitigation strategies.

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree.

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GCCV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GCCV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GCCV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Global variation in diversification rate and species richness are unlinked in plants.

Species richness varies immensely around the world. Variation in the rate of diversification (speciation minus extinction) is often hypothesized to explain this pattern, while alternative explanations invoke time or ecological carrying capacities as drivers. Focusing on seed plants, the world's most important engineers of terrestrial ecosystems, we investigated the role of diversification rate as a link between the environment and global species richness patterns. Applying structural equation modeling to a comprehensive distribution dataset and phylogenetic tree covering all circa 332,000 seed plant species and 99.9% of the world's terrestrial surface (excluding Antarctica), we test five broad hypotheses postulating that diversification serves as a mechanistic link between species richness and climate, climatic stability, seasonality, environmental heterogeneity, or the distribution of biomes. Our results show that the global patterns of species richness and diversification rate are entirely independent. Diversification rates were not highest in warm and wet climates, running counter to the Metabolic Theory of Ecology, one of the dominant explanations for global gradients in species richness. Instead, diversification rates were highest in edaphically diverse, dry areas that have experienced climate change during the Neogene. Meanwhile, we confirmed climate and environmental heterogeneity as the main drivers of species richness, but these effects did not involve diversification rates as a mechanistic link, calling for alternative explanations. We conclude that high species richness is likely driven by the antiquity of wet tropical areas (supporting the "tropical conservatism hypothesis") or the high ecological carrying capacity of warm, wet, and/or environmentally heterogeneous environments.

The Implications of Incongruence between Gene Tree and Species Tree Topologies for Divergence Time Estimation.

Phylogenetic analyses are increasingly being performed with data sets that incorporate hundreds of loci. Due to incomplete lineage sorting, hybridization, and horizontal gene transfer, the gene trees for these loci may often have topologies that differ from each other and from the species tree. The effect of these topological incongruences on divergence time estimation has not been fully investigated. Using a series of simulation experiments and empirical analyses, we demonstrate that when topological incongruence between gene trees and the species tree is not accounted for, the temporal duration of branches in regions of the species tree that are affected by incongruence is underestimated, whilst the duration of other branches is considerably overestimated. This effect becomes more pronounced with higher levels of topological incongruence. We show that this pattern results from the erroneous estimation of the number of substitutions along branches in the species tree, although the effect is modulated by the assumptions inherent to divergence time estimation, such as those relating to the fossil record or among-branch-substitution-rate variation. By only analyzing loci with gene trees that are topologically congruent with the species tree, or only taking into account the branches from each gene tree that are topologically congruent with the species tree, we demonstrate that the effects of topological incongruence can be ameliorated. Nonetheless, even when topologically congruent gene trees or topologically congruent branches are selected, error in divergence time estimates remains. This stems from temporal incongruences between divergence times in species trees and divergence times in gene trees, and more importantly, the difficulty of incorporating necessary assumptions for divergence time estimation. [Divergence time estimation; gene trees; species tree; topological incongruence.].

Concordance-Based Approaches for the Inference of Relationships and Molecular Rates with Phylogenomic Data Sets.

Gene tree conflict is common and finding methods to analyze and alleviate the negative effects that conflict has on species tree analysis is a crucial part of phylogenomics. This study aims to expand the discussion of inferring species trees and molecular branch lengths when conflict is present. Conflict is typically examined in two ways: inferring its prevalence and inferring the influence of the individual genes (how strongly one gene supports any given topology compared to an alternative topology). Here, we examine a procedure for incorporating both conflict and the influence of genes in order to infer evolutionary relationships. All supported relationships in the gene trees are analyzed and the likelihood of the genes constrained to these relationships is summed to provide a likelihood for the relationship. Consensus tree assembly is conducted based on the sum of likelihoods for a given relationship and choosing relationships based on the most likely relationship assuming it does not conflict with a relationship that has a higher likelihood score. If it is not possible for all most likely relationships to be combined into a single bifurcating tree then multiple trees are produced and a consensus tree with a polytomy is created. This procedure allows for more influential genes to have a greater influence on an inferred relationship, does not assume conflict has arisen from any one source and does not force the data set to produce a single bifurcating tree. Using this approach, on three empirical data sets, we examine and discuss the relationship between influence and prevalence of gene tree conflict. We find that in one of the data sets, assembling a bifurcating consensus tree solely composed of the most likely relationships is impossible. To account for conflict in molecular rate analysis we also introduce a concordance-based approach to the summary and estimation of branch lengths suitable for downstream comparative analyses. We demonstrate through simulation that even under high levels of stochastic conflict, the mean and median of the concordant rates recapitulate the true molecular rate better than using a supermatrix approach. Using a large phylogenomic data set, we examine rate heterogeneity across concordant genes with a focus on the branch subtending crown angiosperms. Notably, we find highly variable rates of evolution along the branch subtending crown angiosperms. The approaches outlined here have several limitations, but they also represent some alternative methods for harnessing the complexity of phylogenomic data sets and enrich our inferences of both species relationships and evolutionary processes.[Branch length estimation; consensus tree; gene tree conflict; gene tree filtering; phylogenetics; phylogenomics.].

Chloranthus genome provides insights into the early diversification of angiosperms.

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.